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ENERGY SCALING LAW FOR A SINGLE DISCLINATION IN A THIN 

ELASTIC SHEET 

HEINER OLBERMANN 


Abstract. We consider a single disclination in a thin elastic sheet of thickness h. We 
prove ansatz-free lower bounds for the free elastic energy in three different settings: First, 
for a geometrically fully non-linear plate model, second, for three-dimensional nonlinear 
elasticity, and third, for the Foppl-von Karman plate theory. The lower bounds in the 
first and third result are optimal in the sense that we find upper bounds that are identical 
to the respective lower bounds in the leading order of h. 

1. Introduction 

1.1. Setup and previous work. Consider the following setup: Take a thin elastic sheet 
in the shape of a disc, and remove a sector from it. Then glue the edges of the cut back 
together. The (stress free) reference metric of the sheet is now that of a singular cone. In 
other words, the reference metric is flat away from the centre, where the Gauss curvature 
has a d-type singularity. In the literature, such a defect of the reference metric at the 
origin is often called a “disclination”. We are interested in upper and lower bounds for the 
free elastic energy of the sheet with a single disclination. Note that we we do not impose 
boundary conditions or other constraints. 

Disclinations play an important role, e.g., in the modeling of cell walls |29l I18j . where 
they can serve as the vertices in topologically spherical shells, that take on the shape of a 
(smoothed) icosahedron. Such a behaviour has been observed, e.g., for virus capsids [3]. 
Recent interest in disclinations in thin sheets comes from the analysis of carbon nanohlms 
|28j and the modeling of biological growth processes [9] . An extensive review of the role 
of disclinations in the mechanics of solids, not limited to the case of thin sheets, can be 
found in [28] . 

Considering a single disclination of a comparable deficit without boundary conditions, 
one expects that energy will be minimized by a shape that looks like a cone away from 
the disclination, and deviates from the cone on a ball around the disclination whose size 
is comparable to the thickness of the sheet. In the literature, energy estimates for a single 
disclination have so far been ansatz based (as in [291 im ESI), or it has been assumed 
that the deformation is radially symmetric. In the latter case, the dimensionality of the 
problem is reduced from two to one, and ODE methods can be used. It has been shown 
in [23j that under this assumption, in the Foppl-von Karman approximation, minimizers 
do indeed show the behavior described above. Additionally, this paper contains a quanti¬ 
tative estimate on how fast minimizers approach the singular cone as the distance r from 
the apex of the cone increases. 
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In the recent work [52] , we suggested an ansatz how to prove lower bounds for the free elas¬ 
tic energy of a single disclination without the assumption of radial symmetry. We viewed 
the elastic sheet in the deformed configuration as an immersed Riemannian manifold, and 
focused our analysis on intrinsically defined objects, such as the metric and Gauss curva¬ 
ture of the manifold. The main idea brought forward in this work was that estimates for 
integrals of Gauss curvature can be obtained by interpolation between the metric and the 
Gauss curvature. Lower bounds for integrals of Gauss curvature can be translated to lower 
bounds on the bending energy by the isoperimetric inequality on the sphere. However, in 
|26j . we had to modify the elastic energy and make simplifying assumptions on the shape 
of the sheet to put this ansatz to work. 

Here, we report on an improvement of the method of proof, that makes the simplify¬ 
ing assumptions unnecessary altogether. One further benefit of our modified approach is 
that it makes the use of the jargon of Riemannian geometry superfluous. 

1.2. Statement of results. Our first result regards a single disclination in a geometrically 
fully nonlinear plate model. Let Bi := {a: G : |x| < 1} be the sheet in the reference 
configuration. The singular cone may be described by the mapping : Hi —?■ 

y^{x) = -x/l — -|- A|x|e3 . 

Here, 0 < A < 1 is the height of the singular cone, and is determined by the deficit of the 
disclination at the origin. The reference metric on Bi is given by 

=Dy^{xf' Dy^{x) 

=(1 — A^)x 0 X + A^x (g) X 

=Id2x2 - A^x-^ (g) x-^ , 

where x = x/|x| and x"*- = (—X 2 ,xi)/|x|. The induced metric of a deformation y G 
W‘^’^{Bi-R^) is 

9y = Dy'^Dy. 

The free elastic energy ^ —)■ M is defined by 

IhAiy)= [ {\9y-9A\^ + h‘^\D^y\^)dC\ (1) 

JBi 

where dC^ denotes 2-dimensional Lebesgue measure. The first result of the present paper 
is 

Theorem 1. There exists a constant C = C'(A) with the following property: 

27rA^/i^ ( I log LI — - log I log h\ — C ) < min Ih A(y) < 27rA^L^ (I log LI -|- C) 

V 2 / ;/eVE2'2(Bi;R3) ’ 

for all small enough h. 

Before we come to the statement of the other theorems, let us try to explain the improve¬ 
ment of our techniques that allows us to drop the additional assnmptions we had to make 
in |52|. It consists in noting that the function ^et is controlled by both the 

membrane and the bending term, in different functions spaces. Indeed, it is not difficult 
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to see that is bounded in by the bending energy. The control exerted by 

the membrane energy is slightly more subtle. We recall the very weak form of the Hessian: 

detD^w = (w iW2),i 2 - ^(k,iP),22 - ^(|'R',2p),ii forrcEC^. (2) 

Hence, 

3 1 1 

^detT)2yi = (y,i • y, 2 ),i 2 - -{\v,i \^),22 - -(|y, 2 p),ii for y E , (3) 

i=l 

and we see that i® controlled by the membrane term in Thus, by 

interpolation, we can get control over integrals of ™ Sobolev 

inequality in lower bounds for these integrals can be translated to lower bounds for 
\D‘^y\ (see Section]^, and hence for the bending energy. 

As a consequence of Theoremwe will be able to prove a scaling law in three-dimensional 
elasticity. Let arccos : [—1,1] —[0, vr] denote the inverse of cos : [0,7r] —>• [—1,1]. Define 

:= = (a;i,a; 2 ) E Hi \ {0} : 0 < arccos |-^ < \/l — A 27 r| . 

Let M_ := {(xi, 0) : xi < 0}, and let (/? : \]R_ —>■ M be the angular coordinate satisfying 

X = |x|(cos(/j(x),sin(^(x)) with values in (—7r,7r). We define the map tA : \ M_ —)• Hi 

by 

(\ \ ^{x) . ip{x) \ 

Now let : Hi^a x [—h/2, h/2] —x [—h/2, h/2] be defined by 

i^^(x',X3) = {iA{x'),xs) . 

Note that (i-A^)”^ is defined almost everywhere on Hi x [—h/2, h/2]. We define the set of 
allowed deformations to be 

A,A := {y E W1’2(Hi,a X [-h/2, h/2]-, R^) 

such that To e W^’2 (Hi x [-h/2, h/2]-,R^)'j . 

Let W E [ 0 , 00 )) be the stored energy density that is 50(3)-invariant, i.e., 

W{RF) = W{F) for all R E SO{3) , F E , 
and such that it vanishes on SO{3) and has quadratic growth, i.e., 

W{F) > Cdist2(H, 50(3)) for all F E , (4) 

and 

W{F) < O 2 dist2(H, 50(3)) for F £ Fq , (5) 

where Tq is some neighborhood of 50(3) in M3 x3_ 

The energy functional is given by ■ A,A 1^; 

EhAiY) = ]-[ W{DY)dx. (6) 

^ JBi^Ax[-h/2,h/2] 


We will prove 
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Theorem 2. There exists a constant C = C{A, W) > 0 such that 

^h‘^\logh\< min Eh^Aiv) < Ch'^\log h\ 

^ y^Ah,A 

for all small enough h. 

The third theorem is the analogous result to Theorem in the Foppl-von Karman setting. 
Here, the free elastic energy of the deflection (u,u) G x is 

given by 



2 sym.Du + Dv ® Dv + x"*" (8) x'^ 


2 


+ h^\D‘^v\‘^dx. 


(7) 


Again, the parameter A determines the size of the deficit of the disclination. For a heuristic 
derivation of from Ih^A see Section 5 Our third result is 


Theorem 3. There exists a constant C = C'(A) > 0 with the following property: 

27rA^/i^ (I log/i| — 21og| log/i| — C) < min If\{u,v) < 27rA‘^h‘^{\ log h\ + C) (8) 

(u,v)eA 


for all small enough h. 


1.3. Scientific context. In the last two and a half decades or so, there has been a lot of 
interest in the energy focusing in thin elastic sheets, predominantly in the physics and en¬ 
gineering community. The structures under consideration are typically sharp folds, conical 
singularities, or a (possibly very complex) combination of both. The first analysis of the 
scaling law for the free elastic energy in the sense of an asymptotic expansion for a thin 
elastic shell displaying sharp folds can be found in |34) . The so-called “minimal ridge” has 
been analyzed in [20], see also uniEiiiiaiiHi. Structures with stress focusing in vertices 
(the so-called d-cones) were first investigated in [2|, see also laEUTj. To the interested 
reader, we recommend the overview article by Witten [33|, where many more references 
can be found. 


In the mathematical literature, the most closely related articles to the present work are 
|S1I1I21|. In it has been proved that the elastic energy per unit thickness of a “single 
fold” scales with h^^^, building on results from m- In [241 (3), the following has been 
proved: Take an elastic sheet in the shape of a disc, and prescribe the deformation of the 
sheet at the center and at the boundary, such that it coincides with a (singular) conical 
deformation there. In this class of deformations, the elastic energy per unit thickness 
scales with h‘^\ log/i|. 

On a technical level, the works IHIEIIS] treat variational problems for the functional 0> 
where the reference metric g^ is replaced by the flat 2-dimensional Euclidean metric, with 
certain Dirichlet boundary conditions. These boundary conditions are of a very special 
type: They are chosen such that there exists a unique Lipschitz continuous isometric im- 
mersior0 that satisfies them, which possesses infinite bending energy. In the sequel, we 
will call this Lipschitz continuous isometric immersion the “singular reference configura¬ 
tion” . Also, we will call boundary conditions “tensile”, if a singular reference configuration 

By a Lipschitz continuous isometric immersion, we mean a Lipschitz map y : Bi ^ M. with the 
property that Dy^Dy equals the reference metric almost everywhere. 
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exists. The idea of proof in (HJ [3l [M] is that deviations from the singular reference config¬ 
uration are penalized by the membrane energy. The optimal balance between membrane 
and bending energy yields the lower bound in the elastic energy. 

In the setting of the single disclination, there are no tensile boundary conditions - there 
are none at all. 

To the best of our knowledge, the method of proof using tensile boundary conditions 
has so far been the only one for proving ansatz-free lower bounds in variational models for 
thin elastic sheets with energy focusing in folds and vertices (not counting the predecessor 
|26j of the current work). We expect that the method of proof we present here will have 
applications in variational models that are sufficiently closely related to the single discli¬ 
nation. An interesting open question is whether or not our techniques can also be applied 
in settings with a flat reference metric. 

1.4. /i-principle and rigidity. We want to put Theorem in yet another mathematical 
context. Recall the results by Nash [53] and Kuiper HBj, stating that any strictly short im¬ 
mersion y : M —)• of an n-dimensional Riemannian manifold (M, g) can be arbitrarily 

well approximated by a isometric immersion y in the norm. I.e., given e > 0 and an 

immersion y G with Dy^Dy < g everywhere, there exists y G C^(M; 

with Dy"^Dy = g and ||y — y\\co < £■ These results are relevant for the present context, 
since the leading term in the elastic energy the membrane term, penalizes deviations 
from isometric immersions. By the Nash-Kuiper Theorem, there exists a huge amount 
of mutually vastly different shapes that have arbitrarily small membrane energy. In the 
present setting, one needs to show that all of these many degrees of freedom carry large 
bending energy. 

Next, recall the well known Weyl problem from differential geometry: Given a Riemannian 
metric g with positive Gauss curvature on the two-sphere 5^, find an isometric immersion 
y : ^ Note that by the Nash-Kuiper result above, there exist many solutions if 

we only require y £ C^. On the other hand, by a result by Pogorelov [27|, the solution 
to the Weyl problem is unique (up to Euclidean motions) in the class of immersions with 
bounded extrinsic curvature. This class is dehned by requiring that the pull-back of the 
volume form on 5^ under the Gauss map defines a signed Radon measure. For smooth 
maps, the total variation of that measure is just the L^-norm of the Gauss curvature. We 
see that control over curvature “eliminates” the Nash-Kuiper constructions. 

On a heuristic level, our approach can be viewed as a quantitative version of Pogorelov’s 
rigidity result. In order to see how, note that the different results can be formulated as a 
statement on the measure 

ky : U I —y I KydAy , 

Ju 

defined on open sets U, where KydAy is a certain curvature form (in Pogorelov’s work 
and |26|, the Gauss curvature form; for the geometrically fully nonlinear model treated 
here, det T>^?/idx; for the Foppl-von Karman case, dei D^v dx, which is the Foppl- 

von Karman version of Gauss curvature). The crucial part of Pogorelov’s result is that 
the measure ky is identical to the pull back of the volume form on S‘^ under the Gauss 
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map if the latter is a Radon measure. Here, we give certain quantitative estimates for the 
difference ky — ko, where ko is the curvature measure associated to the reference configu¬ 
ration. These estimates allow for the construction of a lower bound of the bending energy 
via a suitable Sobolev type inequality. 

By imposing tensile boundary conditions as in 0 s 121], one eliminates strictly short 
maps from the set of allowed configurations. In this way, the problem described above is 
circumvented. Here, we describe for the hrst time a setting where the constructions in 
the style of Nash-Kuiper are allowed by the boundary data, but shown to be energetically 
unfavorable. 


1.5. Notation and plan of the paper. The canonical basis vectors of will be denoted 
by ei,..., e^. For a vector x G W^, we denote its components by xi,..., For a vector 
V = {vi,V 2 ) G let u-*- := (—U 2 ,ui). Open balls in are denoted by Br := {x G : 
|x| < r}. The set of linear isometries from to is denoted by 0(2, 3). 

For a function / on M^, we use the notation /^j = dx^f for the partial derivatives. 

The volume of the unit ball in M" will be denoted by a;„ := 7r"'/^/F(n/2-|-l). The symbol 
denotes n-dimensional Lebesgue measure, and T-L^ is the A;-dimensional Hausdorff measure. 
The symbol C denotes numerical constants that only depend on A. A statement such as 
“/ < O^” is to be understood as “there exists a constant C = C'(A) > 0 such that 
/ < O^”. The value of C may vary from one line to the next. By g{h) = o{f{h)) and 
g{h) = 0{f{h)), we mean that 


respectively. 


lim sup 

h-5-O 


\9{h)\ 

\fih)\ 


= 0 , 


\m\ ^ 

lim sup - < oo 

h^o \f{h)\ 


This paper is structured as follows: In Section we state a Sobolev type inequality 
for the Brouwer degree, which will help us to translate control over ky into a lower bound 
for the bending energy later on. Theorems 00 and are proved in Sections [3] m and [5| 
respectively. 


2. A Sobolev type inequality for the Brouwer degree 

In the proof of Theorems and below, we will use the Sobolev inequality to get a lower 
bound on the bending term. From smallness of the membrane and bending energies, we 
will have for most r G [0,1], a lower bound on in the geometrically fully 

nonlinear model, and on 'H‘^{Dv{Br)) in the Foppl-von Karman model (see Propositions 
[^and|^. The Sobolev inequality will give us then a lower bound for \D‘^y\ (in the hrst 
case) and \D‘^v\ (in the second) on dBr- 

In the proof, we are going to need some standard facts concerning the Brouwer degree. 
For a more thorough exposition including proofs of the claims made here, see e.g. [Ill- 
Let U C he bounded, / G and y G M” \ f{dU). Let Ayj denote the 

connected component of M” \ f{dU) that contains y, and let p be a smooth n-form on M"- 
with support in Ayj and g = 1. Then the Brouwer degree deg(y, U, f) is dehned by 

= [ f*g, 

Ju 
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where f* denotes the pull-back under /. It can be shown that this definition is independent 
of the choice of /j., and that deg{-,U, f) is integer-valued, and constant on connected 
components of M”\/(917). By approximation with smooth functions, it can be shown that 
the Brouwer degree has a well defined meaning also for / G (7*^(1/; M") and y G W^\f[dU). 
From now on, let us suppose that U C M"" is bounded and open with Lipschitz boundary. 
For / G C^(17;M"'), and an re-form /x with coefficients in we have 


[ deg(y, [/,/)/x = f />. 

Jr" Ju 

This follows immediately from the defining relation for the Brouwer degree by approxi¬ 
mation with smooth functions. If /x is an exact form, i.e., /x = dw for some re — 1 form u 
with coefficients in then 


[ deg(-,[/,/)/x= [ ridcj) 
JR" Ju 

= [ difu) 
Ju 



We also recall that the space BV{W^) is defined as the set oi g £ 
total variation |T)( 7 |(M"') is finite, 

|T» 5 |(M”) =sup j [ gdivhdx:h£C\R'^-,R^),\h{x)\<lfoT 
iJR" 


(9) 


such that the 

allxGM'^l. (10) 


Let denote the set of re — 1 forms on R"- with coefficients in 

A”"^) =1 ^ u}i^...ir.-idxi, A • • • A dxi^_, : 

Wii-in-i e for all h,.. .,in-i G {!,■ ■ ■ ,n}| . 

Note that the space of re — 1 forms at a given point x is canonically isomorphic to via 


dx^(i) A • • • A dx^(„_i) sgn cr e^(„) , 

for any permutation a : {!,..., re} —>• {!,..., re}. When we write |a;(x)|, we mean the 
norm of the re — 1 form uj{x) under that isomorphism. 

By Q and (10), we have 


|L>(deg(-,[/,/))|(M'^) =sup<^ / deg(-,[/,/)divMx :/i G |/x| < 1 

[jR" 

= sup\f /*a; : a; G A”"^), |a;| < 1 


IdU 


\.JdU 

\Df^^\dn^-^ 


( 11 ) 
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where Df^^{x) denotes the restriction of the gradient of Df to the space tangent to dU 
(which exists almost everywhere by assumption). Finally, recall the Sobolev in¬ 

equality for BV functions on M” (see e.g. [22j . Theorem 2.1), 

\ {n-l)/n 

|/r/("-i)dxJ <C{n)\Df\(R^), (12) 

where C(n) = is the isoperimetric constant in W^. 



Lemma 1. For v G C‘^{Bi) and 0 < r < 1, 


f \D‘^v\d'H^ > ( dvr 

/ det D'^vdx 

JdBr V 

J Bf 


1/2 


Proof. We write w := Dv : ili —t- M^. By approximation with smooth functions, we may 
assume w G C°°{Bi;R‘^). By ( [IT| ), we have 


\D{degi-,Br,w))m^)< iDwldn^ 

JdBr 


By (12), we have 


\D{deg{-,Br,w))\{R‘^) > dvr / \deg{y, Br,w)dy\ 


deg{y,Br,w)dy 
det Dwdx 


> I dvr 


1/2 


1/2 


= I dvr 


' Br 


(13) 


(Id) 


where we used the fact deg(a:, Br, w) G 'L, and hence | deg(x, Br,w)\‘^ > \ deg(x, Br,w)\ for 
almost every x, in the second inequality. Combining (13) and (Id) yields the claim of the 
lemma. □ 


3. Proof of Theorem [T] 

Lemma 2. We have 

inf Ih A(y) < 2vrA2/i2 (I log h\ + C) . 

y£W^,2{Bi-,R3) ’ll 1' 

Proof. Recall the definition of y^, 

y^{x) = -x/l — -|- A|x|e3 . 

Let rj G (^“([O, oo)) with 7/(t) = 0 for t < i, r]{t) = 1 for t > 1 and |r/'| < d, |r/"| < 8. We 
set 

yo{x) = r]{\x\/h)y^{x). 
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We compute 


Dyo{x) (g) ^ r]{x/h)y/l — 


1 0 0 


+ Aea (g) X , 


^2 / N ri (x/h) A „ ^ rj [x h) 

D^yo{x) = —— -y^ ® x ® x -;- 

n 


0 1 0 

\/l — A‘^(2x 0x0x + x'^0x0x'^ 

1 


+ (g) X"^ (g) x) + 2Ae3 (g) X (g) X + ]—(g) X"^ (g) X' 

|x| 

rj{xlh) . , I „ I 

+ , ^ Ae3 ®x^®x^ . 

X 


On Bi \ Bh, we have yo = y^, r]'{\ ■ \/h) = 0, and hence Dy^Dy^ = (Dy^)^Dy^ and 
D^yo = |x|“^Ae 3 (g) x"*- (g) x"*-. On Bh, we see from the formulas above that |I?yo| < C and 
\D^yo\ < Ch~^. (To see the latter, one uses |?/'^| < h on Bh-) Hence the energy can be 
estimated by 


hiuo) < [ {C + h^{Ch ^)^))dx + /i^ f |x| ^A^dx 
<27rA^/i^(| log h\+C). 

This proves the lemma. 


□ 


For the following definition, recall the formulation of the Hessian in its very weak form, 
equation (§• 

Definition 1. Let i G {1,2,3}. By del D'^y^, we denote the distribution 

-7, I A {yf)idyi,2 - yf)idyf)2) • 

^ JBi 

Lemma 3. The distributions det D‘^y^, i = 1,2,3, can be extended to Radon measures. 
Denoting the latter with the same symbol, we have 

det D‘^y^ =0 for i = 1,2 
det D‘^y^ =A^7r(5o , 

where So is the Dirac distribution 6o : (p ip{0). Furthermore, we have 


det D'^yf' 


ivti ■ yt2)d2 - 




in the sense of distributions. 
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Proof. For i = 1,2, the form yiidy ^2 ~ yiidyi 2 vanishes almost everywhere. Furthermore, 
we compute 


2/3^1 - yt,idy% =A^ 


xi ( X 1 X 2 , , \ 

TT nrdxi + -—^dx2 

\x\ \ |x|'^ |x|'’ J 


X2 ( X 1 X 2 , , x| 

-dx 2 + -j—^dxi 


X 


X 


X 


(16) 


=A^ 


— ^dx2 — \ — vT^dxi 
\xr \xr 


It is well known that the differential of the right hand side is just — A^27r(5odxi Adx 2 . This 
proves the first claim of the lemma. The claim (15) follows from a simple integration by 
parts, 

^IJb - i\yt2\^')d) + {{ytyt2)d - i\yt\^),2)) dx 

= {detD‘^y^,(f) . 

This completes the proof of the lemma. □ 

For a given deformation y G let Ky be dehned by 


lip) = (^^det D‘^yi,XBpj 

r 3 

= / det D'^yi dx 
3 1 /. 

= yi,l^yi ,2 - 2 /i,ld 2 /i ,2 • 


i=l 


Proposition 1. Let y G with 

\9y - 9 a\^ + h^\D^yfdC^ < Ch ^\loghi . 




Then 


\Ky — vrA^||£i(g;j^^) < log/i|^/^ for all R G [2h, 1 — 5/i] . 


(17) 


Remark 1. IFe could prove Theorem via Propositions^ and also if we made the 
stronger assumption Ih,A{y) < log/i| in Proposition^ However, the weaker assump¬ 

tion will be convenient in the proof of the result on three-dimensional elasticity, Theorem 


Proof. Let ho = ho{y) G [5h, 6h] be chosen such that 


\9y - 9 A\'^dn^ = I \Dy^ Dy - (Dy^)'^Dy^l'^dH^ < Ch\ logh\ . 


IdBh 


'dBh 
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For r G [/lo, 1], we set 


Fi(r) = / Hy{p)dp, 

J ho 

F 2 {r) ={r - /io)vrA^ , 
F =Fi — F2 . 


By these definitions, 


F'ir) = Y] (det D'^pi - det dp. 

Jho 

Now we set 

Uij = (Dy^Dy - {Dy^f Dy%y . 

By the assumption y G and Lemma we have 

1 1 


F{r) = dp I 

J ^p 


'^ 12,12 — 2 ^ 11,22 — 2 ^ 22,11 


d£2. 


Let k,l G {1,2}, and / G {un, U 12 , U 22 }- We repeatedly apply Gauss’ Theorem, 


I J dp/ ^f^kixjdn^ix) 


'ho JdBp \X 
XI 


Br\BhQ 


\X\ 


/,fc(2:)da; 


(18) 


Xk xi 




\X\ X 


f{x)dn^{x) - 


Br\BhQ 


XI 


f{x)dx 


,k 


We will estimate the norm of the function G : r ^ dp f^kiix)dx. By the triangle 
inequality and (18), we have 


I|G'||li(/io,r) ^ 


r-R 


'ho 


+ 


+ 


r 

'ho 

r 

'ho 


/ ^^fix)dnHx) 

IdBr FI FI 

Xk Xl 


dr 


/ ^^f(x)dn\x) 

IdBf,,^ \x\ kl 

wt) f{x)dx 

lBr\Bh^ VFI/,fc 


dr 


dr 


We estimate the terms appearing on the right hand side using the Cauchy-Schwarz in¬ 
equality. For the first one, we have 

rR 


'ho 


/ ■^:^f{x)dn\x) 

IdBr fI fI 


dr < / 

Xk Xl 

J Bji\Bhg 

F F 


d£"^(x) 


1/2 


< 


F Bjt\Bl^g 

11 
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Secondly, we have 


r 

[ pL^/(x)d7/'(x) 

dr = (i? — ho) 

[ ^^f{x)dn\x) 

Iho 

JdBho 1^1 FI 


JdBho 1^1 


1/2 


<{R-ho)[ \f\^dn^] (271/10)'/' 


IdB, 


hQ 


For the third term, 

rR 


'ho 




Pr ) f{x)dx 


dr <{R — ho) 


I/(®)L^2 




d£"^(x) 


<{R-ho)([ \f\^dcA ([ |x|-2d£2(x)^ 

\Jbu\b^^ J \Jbr\b^^ J 

( \ 1/2 

{R-ho)i[ l/l'd/:^) |log/io|'/' 

\JBR\Bh^ J 


< R- 


Hence, writing |u| := \Dy^ Dy — (Dy^)"^Dy^\, we get 


r-i? 


'ho 


\F{p)\dp<iR-ho) 


[ f \u\‘^dR^ J (27r/io)'/' 

\JdB,^ J 


+ 


+ 


lul^dC^ 


' Rr\R ho 


1/2 


1/2 


1/2-1 


\x\-‘^dC^ 


' Br\R ho 


Br\B ho 
<ChR\ log 

Furthermore, we have the estimate 

rR 


\u\‘^dC^ {tt{R^ - /i^))i/2 


'ho 


\F'\p)\dp = 


' BnXBho 


^ det D'^yi 


i=l 


dx 


< 


D^y\^dx 


J Bi-^h\Bsh 

<C\logh\. 

Now the claim of the proposition follows from the standard interpolation inequality 

see e.g. m, Theorem 7.28. This completes the proof. □ 

Proposition 2. Let yi G C‘^{Bi) for i = 1,2,3. Assume that a > 0 and 

\\Ky — TrA^ll^qg/j < C/i^/^i?^/^| log /i|" for all R G [6/i, 1 — 5/i]. 
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Then 


> 27rA^ (I log/i| — 2alog I log/i| — C) 

j=l -'^l-ShV-Bsh 


Proof. We set 

and choose J G N such that 


ho = 5h\ log h\ 


2'^ho < 1 - 5/i < 2''+"ho . 


2a 


J +1 2 


Note that this choice implies Jlog2 —C < | log ho\ < (J+1) log 2+ (7. For j = 0,..., J+1, 
we will write Rj := 2^ho- 

We use Jensen’s inequality to get for r G [5/i, 1 — 5/i], 

^ -.2„. |2^ n_,. ^dBr 


IdBr 


\D"‘yifdTL^{x) > 27rr 


2'Kr 


for i G {1, 2, 3}. 


Recall the definition of Ky in ( |17[ ) . We use Lemma on the right hand side in (20), 
over i, integrate over r, and use the triangle inequality to obtain 


( 20 ) 


sum 


E 


^l—bh 




D yi\ dx>— I —4Tr\Ky{r)\ 


' bh 


Reducing the domain of integration and using again the triangle inequality, we have 




iDW^dx >27rA2 ^ - 2 T" -\Ky{r) - ttA^ 

J Ro ^ J Ro ^ 

=27rA^ (I log h\ — 2a log | log h\ — C) 

-2/" -ttA^i . 

jRo ^ 


( 21 ) 


Using the assumption (19), we may estimate the last term on the right hand side as an 
error term, 


'■R.j 


' Ro 


dr, 


J f-R 


■3 (It- 


— |K,,(r) — ttA^ 


j = l '' Bj-1 

J 

- X/ 11 LI (/i,Hi) 


i=i 

i=i 


where in the second inequality, we have used (46), and the last inequality is just the 
summation of a geometric series. Inserting this in (21), the proof of the proposition is 
complete. □ 
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Proof of Theorem\^ First note that the minimum is actually attained by some y G M^), 

since the functional 0 is coercive and lower semi-continuous. The upper bound holds 
by Lemma[^ Assume that y G M^) with Ih,Aiy) < 27rA^/i^ (| log /i| +C). By 

density of in we may assume y G By Propositionthe assumption 

of Proposition 2 holds with a = |. The lower bound now follows from that of Proposition 

H □ 


4. Three-dimensional elasticity - Proof of Theorem [2] 

This section is concerned with the transition from two dimensions to three. The three- 
dimensional model is defined by the free elastic energy functional (§. We first consider a 
special class of deformations. Let y G 1^^) be an immersion, and let the three- 

dimensional deformation Y be given by the following Kirchhoff-Love ansatz: 

Y{x', xs) = y{x') + xsVy for {x , X 3 ) G Bi^a x [-h/2, h/2] , 

where Uy = y^i A y^ 2 /\y,i A y^ 2 \ denotes the unit normal. With the assumption ([^ on 
the stored energy density W, and provided that ||Z)y||i;,oo < C, it is not difficult to show 
(cf. the proof of Lemma below) that 

/ W{DY{x', X3))dx3 < C (dist2(Dy(x'), 0(2, 3)) + h^\D^y{x^)\^) , 

J [-h/2,h/2] 

and hence 


Eh,A{Y)<C dist\Dy,0{2,3)) + h^\D^y\^d£\ (22) 

J 

Whenever a three-dimensional model and a plate model are related to each other through 
this estimate for the Kirchhoff-Love ansatz, we say for short that they correspond to each 
other. 

As a matter of fact, it is also possible to obtain lower bounds for three-dimensional mod¬ 
els from lower bounds in the corresponding plate models. This is based on the Geometric 
Rigidity Theorem by Friesecke, James, Muller, which we cite now: 

Theorem 4. [13] Let n > 2 , and let U C M"" be a bounded connected Lipschitz do¬ 
main. There exists a constant C* = C*{U) with the following property: For every 
u G there exists R G SO{n) such that 

\\Du — R\\i, 2 ijj-^ < C* f dist"^{Du, SO{n))dx . 

J TJ 

The constant C* is invariant under rescalings ofU. 

The deduction of lower bounds in three-dimensional elasticity starting from lower bounds 
for the corresponding plate models via Theorem]^ has been carried out in [30] . The main 
strategy of our proof is going to be similar to the one given in that reference. Since the 
estimates needed here differ in many details from those in [30] . we nevertheless include a 
full proof here for the convenience of the reader. 

The idea is as follows: From a given deformation in the three-dimensional model, one 
constructs a deformation in the two-dimensional one through averaging and smoothing. 
On small boxes of sidelength comparable to h in the 3-dimensional domain, the defor¬ 
mation gradient is close to 50(3) by Theorem This implies that the gradient of the 

14 
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2-dimensional deformation derived from it is close to 0(2,3). In this way, the membrane 
term of the two-dimensional model can be estimated from above by the elastic energy in 
the 3-dimensional one. The bending term can be dealt with in a similar way. 

We now recall some notation, and introduce some more. 


Recall from Section 1.2 that i/s. : Hi ,a Hi is defined by 

ip{x 


i/s{x) = \x\ cos 


, cos 


ip{x) 


Vl - A2 ’ Vl - A2, 
where <^(x) is the angular coordinate of x. To alleviate the notation, we write 

i(x) = iA(x) 

for X G Hi,A- When convenient, we will use the notation z := r(a:). 

On r(Hi,A) = Hi \ {(xi,0) : xi < 0}, i has a well defined inverse, that we denote by 

j : Hi \ {(xi, 0) : xi < 0} Hi ,a . 

For y G let y : A/i — )■ be dehned hy y = y o j. 

In some places, it will be convenient to work, instead of Hi, a, with the smaller domain 


Wh,A := S X G Hi_ 5 fe \ Hs/i : 0 < arccos ^ < \/1 - A'^tt 


x| 

By slight abuse of notation, we will not distinguish between the map i, dehned on Hi, a, and 
its restriction to In the same way, we will not distinguish between j : Hi \M_ —)• Hi,a 

and its restriction to := i{u}h,A) = \ H 5 / 1 ) \ {(xi, 0 ) : xi < 0 }. 

If for two quantities f,g, there exists a constant C = C'(A) such that 

y<g<Cf, 

we will write f — g. 

Furthermore, for tensors a G m 3 >< 2 x 2 ^ ^ m 2 x 2 ^ ^ g ]^ 3 x 2 x 2 dehned by 


(ci . ^ ^ Oiilm^ljbmk 

/,m=l 

for i G {1,2,3}, j,k G {1,2}. 

4.1. The corresponding plate theory. To prove Theorem we will hrst have to prove 
an energy scaling law for a two-dimensional plate model that slightly differs from Q, 
namely, the one on the right hand side in (22). The reason why we do not choose to work 
with the three-dimensional model corresponding to Q is that we want to apply Theorem 
and the latter requires a hat reference metric - not only in the three-dimensional model, 
but in the corresponding two-dimensional one too. 


:= 


■= {y e : y o i"' G . 
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Our two-dimensional plate model is given by Jh,A ■ —>• 

JhAiy) = [ \9y- Id2x2|^ + . 


'Si 


It will be convenient to define a slight variation of Jh,Aj that consists of a reduction of the 
domain of integration for the energy density. Let 

JhAiy) = [ \gy -Id2x2\^ + h^\D‘^y\‘^dC^ . 

Lemma 4. With y := y o j, we have for every y G 
JhAiy)- \gy{z) - gAiz)f 


+ h^ 


D^y{z) —^ {Dy{z) ■ z) ® z^ ® z-^ d£^(z), 


(23) 


JhAiy)- \gyiz) - gAiz)\" 

J Bi 


+ h^ 


D^y{z) - iDy{z) ■ z) ® z^ ® z^ d£^{z). 


Proof. For x G cohA^ compute 


Di\x =z (8) X -h 


v/1-A2 


z (8) X 


D^i\x = - 


A2 


(24) 


|x| (8) x"*" ® X"' 


1-A2' 

The hrst equation implies 

Dl\x zz Id2x2 (25) 

in the sense of positive dehnite matrices. By the chain rule, Dy\x = Dy\^^x) ' Di'\x for 
y G A"^^ and x G iOhA ■ a consequence of the second equation in (24), we have 

D y\x =D y\i^[x) • iJ^d\xi J9i\x) + J^y\i{x) ' J^ '-U 


A 

=D‘^y\i{x) ■ iDi\x, Dl\x) - (-C>i/L(x) • 5) <8) x-^ 0 x^ 


Hence, 


JhAiy) = I DL\lDy\'f^^.^Dy\,(^x)DL\x - Id2x2 


'^h,A 




^ y\i,(x) ■ iJ^l'lx: J^l'lx) 


(26) 


A2 


\x\ ^ iDy\,(^x)-z) ®X^ ®X^ 


1 - A2' 

By (25), we have for every a G and for every b G 

|o : (L>r|“\L>i|“^)p ~|o|2 

iDi\f^)'^ zz\h\\ 
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Using this in (26), we get 




JhAiy) - / ^ (^4 


‘^h,A 
2 


+ h 


D‘^y\i(x) - (l)yL(x) ■z)®z^®z^ dC‘^{x) 


(27) 


Multiplying this relation with a factor |det Dl\x\ = (1 — ~ 1, and using the coarea 

formula, we obtain the first relation in ( |23[ ); the second is obtained in exactly the same 
way. □ 

Lemma 5. There exists G with 

\Dy^\ <C 
\DuyH \ <C|Z)Vl 
JhAiy"") <c/i2|iog/i|. 

Proof. Let : i?i —>■ be the test function that we constructed in Lemma i.e., for 

X G Bi, 


y^(x) = y{\x\/h)y‘^{x ), 


where rj G C°°([0,oo)) with r]{t) = 0 for t < ^, r]{t) = 1 for t > 1, |y'| < 4, \r]"\ < 8; and 
is given by 

y^{x) = -x/l — A^x + Aesl^l. 


We set y^ := y^ o i. From the explicit formulas for Dy^, D'^y^ in the proof of Lemma 0 
and (24), we see that \Dy^\ < C, Dy^ = 0 on \Duyh\ < C\D'^y^{z)\ < C\z\^ 

for z G Bi\Bf^/ 2 - Furthermore, we have gyh = yA on Bi \ B^ and \gyh \ < C on Bfi. Hence, 
using (23), we may estimate 


JhAiy’^) - / \9yh{z) - gA{z)\^ + 


'Bi 


D'^y'^{z) — ^{Dy'^(z) ■ z)'S> z'^'S> z-'- d£^(2;) 


<[ CdC‘^ + h^[ C\z\-‘^dC^{z) 

<C/l2(|log/l| + 1) . 


This proves the lemma. □ 

Proposition 3. For all h small enough, we have 

]^hf‘\\ogh\< inf JhAiy)- 
17 
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Proof. Assume that Jh,A{y) < Ch'^\logh\ (otherwise there is nothing to show). Setting 
y := y o j, we may estimate 




z\-^\Dy{z) ■ z\^dC^iz) 

< [ \z\-‘^dC^{z) + 

<27r(| log/i| — C) 


+ 


Bi-5h\B^h 


1/2 


[\Dy{z) ■ z\'^ - l) dC^{z) 


1/2 


(28) 


\gy - gA?dC^{z) 


, , kr"d/:2(z) / 

\JBi_^i^\B^f^ J Si_5(,\S5h 

<27r| log/i| + C [h^\ log/i|)^^^ 

<27r| log/i| + o(| log/i|), 

where we have used the Cauchy-Schwarz inequality and Lemmain the second inequality. 
Again by Lemma and Young’s inequality, we have 


' Bi-5h\B5h 


\D^y\dC‘^ <Ch^ ( I log/i| + [ \z\-^\Dy ■ z\^dC^{z) 

\ J Bi-5h\Bsh ^ 


<Ch^\\ogh\. 


Hence, by Propositionthe assumption of Propositionis fulfilled for y with a = |. By 
Proposition 


' BishXBsh 


|Z)i/P > 27rA^ ( I log/i| — - log | log/i| — C 


By (28), (29) and the triangle inequality, we have 


D^y—lA^\z\ ^ {Dy ■ z) ® z-^ ® z-^ 


L^{Br.5h\B5h) 


> (^\/27rA2| log/i| - v^27rA4| log/i|^ -o(|log/i|) 
=27rA(l — A)| log/i| — o(| log/i|). 

Combining this last estimate with Lemma proves the proposition. 


(29) 


□ 


4.2. Reduction of domain and covers by boxes. As we explained at the beginning 
of the present section, the core step in the passage from two to three dimensions is to 
cover the three-dimensional domain by small boxes, and to approximate the deformation 
by a Euclidean motion on each of these boxes. In this approach, one needs to take special 
care of the boundary of the two-dimensional domain. Part of that boundary (namely, 
dBi /\ n dBi) has already been taken care of by “shrinking” i?i,A to uJh,A in the two- 
dimensional plate theory, as we did in the previous subsection. The present subsection’s 
aim is to control the boxes that are close to the remaining part of the boundary. The 
proof of this subsection’s main statement (Lemma below) is straightforward; however, 
we will need to introduce a certain amount of auxiliary notation. 
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We will be interested in coverings of uJh^A by small squares of sidelength h or 3h. For 
a G and r > 0, let Qa,r '■= a + [—t/2, r/2]^. We define 


Qh :={a G hZ^ : Qa,h n u}h,A / 0} 
Qh • = {“ ^ Qh ■ Qa,3h ^ ^1,a} ' 


(30) 


Furthermore, let 

-Bi,A := 


:= |x G \ {0} : (l - \/l - A 2 ) 


Xl 

TT < arccos ^—r < TT 

X 


(±) . 


^ • -l^i,a ^ -Bi,a that translate the cut in the domain in 


B, 


Now we will define maps 
angular direction. Let (pA '■= (1 — Vl — A2)27r and let the rotation Sa G 50(2) be defined 
by 

cos Pa — sin pA 


Sa = 


sin pA cos Pa 


Let 


V := cos 


1 \/l — A^tt, ± sin Vl - AStt) . 


Note that are chosen such that [0, u+)U[0, v~) = dBi^A\dBi, where by [a, 6), we denote 
the line segment connecting a with 6 in including a, but excluding b. Furthermore, 
note that Sav~ = u"''. 

Now let 

•Bi,a,+ ■={^ ^ -®i,A : X ■ (u+)-^ > 0} 

•Bra■={^ ^ Bi^a : x ■ { v ~)-^ < 0} 


and define /a,± by 


if X G Bi_a,+ 


fA,+ ix) = 


fA-ix) = 


Sax if X G Bi^a,- 


X 


else 


For a sketch of Bi^a,± and /a,±, see Figure 
Note that for x G Bi a, 


S^^x if 5^^x G Bi^a 


- 1 , 


/a,- ° /a,+ (x) = 


S^^x else. 


For Y G Ah,A, we define y(=*=) : .Bi_a x [—h/2, h/2] —)> by 


(31) 


y(=^)(x',X3) = y(/A,±(a:')>a:3) ■ (32) 

The main statement of this subsection is 
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Figure 1. The maps /a,±- Colored in gray; in the top right, the set 
; in the top right, /a,+ (-Bi,a,+ ); in the bottom left, i?i^A,-; in the 
bottom right, /a,-(-Bi^A,-)- 


Lemma 6. Let Y G ^h,A- Then 


V / dist2(Dy(±),50(3))d£3 

<3/" dist2(Dy,50(3))d£^ 

■tBi,Ax[-V2,V2] 


(33) 


Proof. First, we note that as a direct consequence of the definition of Q*^, we have that 
for every a G Qa,?,h C Bi^^. Next, for g G L^(5i,a), let 


5'A,± =50 /a,± • 

Note that away from their discontinuity sets (which are null sets), /a.g are isometries. 
Hence, we have 



|5A,±|d£2 <3 f 

J Li,A 


\9A,±\dC^ 


=3 



(34) 
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For every x' G i?i,A that is not on the discontinuity set of /a,±) there exists R G SO{3) 
such that we have 

DY^^Hx^,X 3) = DY(fA,±{^'),X3)R, 

and hence for every (x',X 3 ) G -Bi^a x [—h/2,h/2], 

dist2(ByW(x',X3),BO(3)) = dist2(By(/A,±(x'),X3),SO(3)). (35) 

For fixed X 3 G [—h/2,h/2], combining (34), (35) and the fact that the discontinuity set of 
/a,± is an null set, yields 

f dist‘^{DY^^\x',X3),SO{3))dC‘^{x') <3 f dist‘^ {DY (x', X3), S O {3))dC‘^ {x'), 

a£Q^ 

where we have used the coarea formula to make the change of variables x' = fA±ix') on 
the right hand side. Integrating in X3 yields the claim of the lemma. □ 


4.3. Proof of Theorem | 
Lemma 7. IFe have 


inf Eh^A{Y) < Ch^\ log/i| . 
y ^Ah,A 


Proof. Let : 0Jh,A 1^^ be as in the statement of Lemma Now, for (x',X 3 ) G 
Bi,a X [—/i/2,h/2], we set 


Y^(x',X 3) = y^{x') +X3V h{x') 


We compute 


DY^{x', X3) = Dy^{x') + i^yh{x') (g) 63 + X3DUyh{x') 


^'^{x\x3) = Dy'^{x' 

By definition of the unit normal Vyh , we have 
1 - 2 / 


d\sY{Dy^{x') + Vyh{x') ® 63, SO{3)) < dist^(-Dy^(x'), 0(2, 3)). 

By Lemma 5 we have \Dvyh\^ < Ol-D^y^p, and hence, we get by integrating in X3, 


-V2,/i/2] 


disF^(-Dy"(x', X3), SO{3))dx3 


< C (dist\Dy^{x'), 0(2, 3)) + h^\D^y^{x')\^) 


(36) 


for all x' G u}h,A- 

Since W{F) < C dist^(F, 50(3)) in a neighborhood of 50(3) by ([^, we have that for any 
compact set F C there exists a constant C = C{W,F) such that 


W{F) <Cdist^{F,SO{3)) for F € F. 


Since \Dy^\ < C, we get 


W{DY^) < C{W) dist2(B/, 50(3)). 
21 
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Hence, using (36), 


EhAiy"") = f dx' / dx3W{DY'^) 

"'-Bi.a J[-h/2,h/2] 

<C{W) f dx'(^dist2(I)?/'*, 0(2,3 ))+ /i 2 |dVP) 


' J^l,A 
,2 


<C{W)h^\logh\ , 

where the last estimate holds by Lemma This proves the lemma. □ 

Proof of Theorem^ The upper bound has been proved in Lemmawhile the existence 
of a minimizer follows from coercivity and lower semicontinuity of in x 

[—h/2, /i/2];]R^). It remains to show the lower bound. 

Step 1. Extension of the domain and mollification. Assume that Y G x 

[—h/2,h/2];M.^). We define —)• by setting 

(L/'Vi ■= 


y(P(x'):=f y(x',X3)dx. 

J [-h/2,h/2] 

For each x' G Bi^^, we set 

y(b±)(x') := y^'^\fA,±ix')) . 

Next, we mollify; let rj G 0“(M^) with fj^ 2 V = 1 and suppr/ C B^. We set %(•) := 
h-^T/(-/h). Now, for x' G UasQ^VQ* 

y(^)(x') = {yh * y^^''){x '). 

For x' G we define 

= [yh * y^^'^'^){xf). 

Since linear transformations and convolution commute, we have by (31), 

= y^‘^’~\S'^^x') for all x' G Hi^a with x' G Hi ,a • 

Now let 


(37) 


Wh,A,± = < X G [J Qa,h n U}h,A ■ fA,±{x') = x' 

I “62); 

Note that we have 

^h,A — ^h,A,+ U ^h,A,— U Qa,h j 

where the sets on the right hand side are mutually disjoint. On uJh^A,+ define 

1 /( 2 ) by 

y(2)(a:;') = ?/(2>±)(x') if x' G 0Jh,A,± ■ 

By (37), we have 

y(^) oj GW^’\Bi_5h\B5h-,M.^). 

Step 2. Application of Theoremon small boxes. Let a G Qh \ Q^- By Theoremthere 
exists Ra,h 2 SO{3) such that 


lQa,3hX[-h/2,h/2] 


dist2(Hy,50(3))dx < C 
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For a G there exist G SO{3) such that 


/ dist2(F)yW,50(3))dx < C / 

lQa,3h>i[-h/2,h/2] J Qa,3h>i[-h/2,h/2] 


|I)y(±)_^(±)|2da;, (39) 


where has been defined in (32). For a G Qh \ Q^, let the 3 by 2 matrix consisting of 
the first two columns of Ra,h be denoted by Ra,h- For a G Qh, let R^^^ denote the first 


two columns of respectively. 

Step 3. Conclusion. Note that for every a G Qh\Q*h, and every x' G Qa, 3 h, 


Dy^^\x')-R, 


’a,h 


< 


-h/2,h/2] 


\DY{x',X 3 ) - Ra,h\ d2;3 ■ 


An analogous formula holds for a G Q\. Now we have 



dist2(F>y(2)^0(2^3))d£2 < 


E 


d^QhXQt 


' Qa 




+ Y. ! dist2(Z)y(2,d^o(2,3))d£2 

aes;|;,ie{±} 


< Y [ \Dy^^^ - Ra,h\^dC^ 

aeQft\S* 

+ E / loy*’- 

aes;;,ie{±} 


(40) 


We use the definition of and Young’s inequality for convolutions to get (for a G Qh\Q% 



with analogous estimates for a G Q\. Inserting in (40), we get 


23 







H. OLBERMANN 




' ‘^h,A 


<c( Y. [ \Dy^^'> - RaM^dC^ 

VaeQh\Qr 

+ E f --Rl‘lPd£" 


aes^,ie{±} ■'*^“’3/1 


<Ch' 




\DY - Ra,h\^dC^ 


V «6S/i\QJ;fe/2,h/2] 


+ V 

aeS^,ie{±} ■^Qo;3hX[-h/2,h/2] 


|E>yW _ |2^^3 


dist2(Z)y,50(3))d£^ 


+ E / 

aeQ^,ie{±} "'Qa,3/iX[-/i/2,/i/2] 

Using Lemma in (41), we get 


aeQhXQl'^ /i/2,/i/2] 

dist2(Dy«,50(3))dE3 ] . 


(41) 


dist2(Zly(2),0(2,3))d£2 < c/i"^ 


' ^h,A 


'Bi_AX[-/i/ 2 ,h/ 2 ] 


dist2(^ly,50(3))d£^ (42) 


We come to the estimate of the bending term, using again the definition of and Young’s 
inequality for convolutions: 


f |Wfd£2< ^ f ^Dyh*{Dyd) 

J^h,A a£Qh\Qt'^^°-’^ 


d £2 


+ E 


Qa,h 




d£^ 


<Ch 


-2 


^ L 

&Qh\Ql 


Dy^^^-Ra,h 


+ E 

aGQ^,iG{=b} 


Qa,3h 


Os“'‘’ - <1 


d£^ 


(43) 


Using the chain of inequalities from (41) starting from the second one, and again Lemma 
we get from (43), 


\D^y^^^\‘^dx < Ch 


-3 


' ‘^h,A 


'Bi AX[-h/2,h/2] 

' 24 


dist2(Dy,50(3))d£^ 


(44) 
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Combining (42) and (44), we get 


JhAiy^^^) < ch-^ 

By (§, this implies 

JhAiy^^^) < Ch-^ 


'Bi.Ax[-h/2,/i/2] 


dist2(Dy,50(3))d/:^ 


W{DY)d/:^ = CEhAiY) 


'-Bi,ax[-/i/2,/i/2] 

By Proposition!^ this proves the lower bound, and completes the proof of the theorem. □ 


5. The Foppl-von Karman setting - Proof of Theorem [s] 

The Foppl-von Karman (FvK) plate model is a very popular model for phenomena that 
include large deflections. Strictly speaking, there exists no rigorous justification for the 
validity of the model when large deflections occur. The FvK model can be derived from 
three-dimensional nonlinear elasticity as a low energy Gamma-limit, see |13( 112) . This 
rigorously justifies the FvK model, but only as a model for small deflections. 
Nevertheless, the FvK equations, first formulated more than hundred years ago [35], have 
a long and quite successful history of describing large deformation phenomena in engineer¬ 
ing, from the design of submarine hulls over paper crumpling and wrinkling 

to the mechanical properties of cell walls [291 ITS] . In the context of stability and buckling, 
several interesting points on why the FvK equations are relevant, even in regimes that 
they seemingly do not describe, have been made in m- 

In the applications mentioned above, the FvK model usually still contains a small pa¬ 
rameter h that is interpreted as the (rescaled) thickness of the sheet, whereas in the 
Gamma-Limit from [12], the thickness has disappeared from the model. Strictly speak¬ 
ing, it is unclear if plate models that still contain the small parameter h have a rigorous 
meaning asymptotically. 

Theorem identifies the energy scaling of a single disclination in the FvK model that 
includes the small parameter h. The size of the disclination is of the same order as the 
“natural” deflection size in the FvK model, and is determined by the parameter A. 

To see the relation between I^a let e be a small parameter, and set A = eA, 

h = eh. Assume further that {u,v) £ x and that the deformation 

y : Bi —?• is given by 

y{x) = X + e^u + eve^ . 


Dy =ei (8) Cl -|- 62 <8) 62 -I- e^Du + 663 (g) Dv 
Dy^’^Dy =Id2x2 + e^{Du + Du^ + Dv (g) Dv) + 0{e^) 

Dy^Dy — =e^ (^Du + DvF + Dv (g) Dv + A^x-*“ (g) , 
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and hence 

hKiy) = l \gy-gK? + h^\Dh?dC^ 

=e^ / \2sYva.Du + Dv®Dv +®x^\^+ h?\D^v\^dC^+ 0{e^). 

Jbi 

Hence, we see that Ih\{u, v) is the lowest order term in e of Ieh,eA{x i—)■ a: + e^u + eve^). 
We remark that the above expansion in e can not be used to deduce Theorem from 
Theorem[^in a trivial way. Nevertheless, the proof of Theoremj^is very similar to the one 
of Theoremj^ The role that was played by the control variable Ky{r) = det D'^yidx 

there will be played by 

K^^(r) := f detD'^vdx (45) 

J Br 

here. We recall that detD'^v is the Foppl-von Karman version of Gauss curvature - so 
our control variable is a curvature integral, and by Gauss’ equation may be thought of as 
the (oriented) volume of balls Bj. under the push-forward of the Gauss map. 

The only significant difference between the two cases is that we will have to identify the 
right terms in the membrane energy to make estimates for in the norm. 


We start off with the upper bound: 

Lemma 8. We have 

inf /)(^(rt, u) < 27rA^/i^l log/i| + C/i^ . 

(u,v)&A ’ 

Proof. Let rj G C°°([0,oo)) with r]{x) = 0 for x < ^, r]{x) = 1 for x > 1 and \r]'\ < 4, 
\r]"\ < 8. Now we define : Hi —)■ : Hi —)• M by 

■= - —g{\x\/h)x 
v^{r,(p) ;=Ar 7 (|x|/h)|x| . 

The membrane energy density |2symHrt^ + Dv^ (8> Dv^ + A^x"*" G x“’'|^ vanishes on Hi \ 
Bh, while it is bounded by a constant C that is independent of h on B^. Hence 

2 

2 sy m D ii!^ + Dv’^ (g) Dv^ + A^x-^ (g) x^ 


Hmembrane • — 

Furthermore, 




' dx < Ch^ . 


= A (^^g''{\x\/h) + ^??'(|x|//i)^ X (g) X ^^7?'(|x|//i) -k r/(|x|//i)^ x-^ (g): 


Hence 


<Ch-^ for |x| < , 

|x| ■ for |x| > h, 


| T ) 2 ^ h |2 ^^ 2|^|-2 


h ^Hbending := [ jH^x^^dx < f Ch ^dx + 27r f A^ — 
J Bi J Bh Jh X 

<27rA^| log/i| + C. 
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This proves the lemma. □ 

We come to the interpolation between metric and curvature that will yield an estimate 
for kI^. 

Proposition 4. Let {u,v) G with < 27rA^/i^(| log/i| + C). Then 

ll<^ - '^^‘^\\Li{h,R) ^ logh\ . (46) 

Proof. In the present proof, we will use polar coordinates r, ip on Bi. The unit vectors in 
r and (p direction are denoted by e,. = x and e^ = x-^ respectively. The vector field u will 
be written as u = UrCr + u^e^p. The membrane energy in these coordinates is given by 


^2-k 


rdr / dipi \2ur^r + + \2r ^{nP^ + Ur) + {r + Li. 


- 1 „. , a 2|2 


'0 ~'o 


T 2 “t“ V Rp T j • 

For r > 0, we have 

I detD‘^vdx=- I vpdvp — vpdvp. 

J Br 2 Jqb^ 

With = v^r cos ip — sin ip and vp = v^r sin ip + r~^v^p cos ip, we compute 


(47) 


vpdvp — vpdvp =adr + + 




dip 


{v,rV,p) 


=adr + {v^rY - — + 




dip, 


(48) 


where a is some function of r, ip that will be irrelevant for our purpose. For h < s,r <1, 
let 

ns n 27 r 

a{s) := dr dip {2ur^r + v^r) 


’ h Jo 

r-2TT 


V , 


b{r) := J dip ^2uj. H —— + 


Now we set F{s) := ^ {a{s) — b{s)), and compute 




Comparing this last expression with (47) and (48), we have 


F'(r) 

F"(r) 


= - vr A^ + 


det D^vdx 


' Br 


= J rdip det D^v . 
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This gives us the needed bound on 


r2-K 


W'\\L\h,i)< rdr difldetD 


'0 


< / \D‘^v\dx 
Jbi 

<C\logh\. 


(50) 


The bound on F works as follows, 






|6(r)|dr = / dr 


d(/j ( 2ur H —— + A^r 

V 


< 


rdrd(p\2r ^+ Ur) + {r 


1/2 / rR \ V2 
rdr 


and 


<C/i|log/i|^/2^, 

\a{s)\ = 

from which we get 

In conclusion, 


dr J d(p (2Ur,r + v^j.) — (^J J |2rtr,r + 

<Ch\ log/i| 


V2 / ..^y/2 
1 '■ / 


r-R 


|a(r)|dr < ChR\ log h\ . 


fR 

/ |F(r)|dr < C/iii|log/r| . 

Jh 

Using the standard interpolation inequality 
and (|49|), we obtain the claim of the proposition. 


(51) 


□ 

The proof of Theorem now almost works exactly in the same way as the proof of Theorem 

m 

Proof of Theorem The existence of a minimum follows from coercivity and lower semi¬ 
continuity of the functional Q. The upper bound follows from Lemma Assume that 
{u,v) G VT^’^(ili; M^) x M^) with < 27rA^/i^(| log/i| +C). By density of 

C'^(i?i;M^) in IT^’^(Bi;M^) x IU^’^(i?i;M^), we may assume {u^v) G C'^(Bi;M^). By 
Proposition!^ we may assume that the assumptions of Proposition!^ hold with a = 1 and 

yi = 1/2 = 0, yz = v . 


With this choice of yi, we have Ky = and 


^ [ \DW^dx = [ \D 

JBi JBi 


^rPdx. 
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Hence, the lower bound follows from the statement of Proposition This completes the 
proof of the theorem. □ 
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